Three-photon states in nonlinear crystal superlattices 
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It has been a longstanding goal in quantum optics to realize controllable sources generating 
joint multiphoton states, particularly, photon triplet with arbitrary spectral characteristics. We 
demonstrate that such sources can be realized via cascaded parametric down-conversion (PDC) in 
superlattice structures of nonlinear and linear segments. We consider scheme that involves two 
parametric processes: ojq — >• uj\ + u)2, ^2 — >■ wi +wi under pulsed pump and investigate spontaneous 
creation of photon triplet as well as generation of high-intensity mode in intracavity three-photon 
splitting. We show preparation of Greenberger-Horne-Zeilinger polarization entangled states in cas- 
caded type-II and type-I PDC in framework of consideration dual-grid structure that involves two 
periodically-poled crystals. We demonstrate the method of compensation of the dispersive effects in 
non-linear segments by appropriately chosen linear dispersive segments of superlattice for prepara- 
tion heralded joint states of two polarized photons. In the case of intracavity three-photon splitting, 
we concentrate on investigation of photon-number distributions, third-order photon-number cor- 
relation function as well as the Wigner functions. These quantities are observed both for short 
interaction time intervals and in over transient regime, when dissipative effects are essential. 

PACS numbers: 42.65.Lm, 42.50.Dv, 42.65.Yj 



I. INTRODUCTION 



Quantum communication and optical quantum com- 
putation rely on controlled preparation of well-defined 
photonic states as well as entangled states of electromag- 
netic modes. Generating of these states is a primary 
task for the application of quantum information process- 
ing. Particularly, the experimental preparation, manip- 
ulation, and detection of multiphoton states are of great 
interest for the implementation of quantum communica- 
tion schemes, quantum gates and for fundamental tests 
of quantum theory [IHH. There exist a number of pro- 
posal for implementing optical quantum-state engineer- 
ing (see, for example, reviewed paper [4] ) . Recently, this 
program has been realized in details for single-photon and 
single-mode states. The method of manipulating overall 
group delay mismatches between interacting in a multi- 
layered structures (superlattice structures) by compen- 
sating dispersive effects in nonlinear segments by chosen 
linear segments have been developed |5|-l7| for synthesis 
of twin photon states. In this way, the demonstration 
of heralded single photons prepared in pure quantum 
states from a spontaneous parametric down-conversion 
(SPDC) has been done in series of papers § and engi- 
neering of quantum-optical state has been demonstrated 
in Q. There are some attempts to expand this strategy 
for producing heralded two-photon polarization entan- 
gled states, since one promising approach to the practi- 
cal realization many of these tasks relies on qubits that 
are encoded in the polarization states of single photons. 



It was shown that the production of one heralded po- 
larization entangled photon pair using only conventional 
down-conversion sources, linear optical elements, and 
projective measurements requires at least three initial 
pairs 10]. An experimental scheme for producing her- 
alded two-photon entanglement, which relies on triple- 
pair emission from a single down-conversion_ source was 
proposed [Til ] and experimentally realized jl~2j ] . 

It has an important goal of quantum optics to re- 
alize effective sources producing three-photon entan- 
gled states. Up to now, several physical systems 
have been proposed for the direct generation of photon 
triplet including cascaded spontaneous parametric down- 
conversion HH, [l4 1 , quantum dots [lj| as well as electron- 
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positron annihilation (l6| . Experimentally, three-photon 
down-conversion was studied in third-order nonlinear me- 
dia [l7l - il9| and also by using cascaded second-order non- 
linear parametric processes [2fj| . Most recently, direct 
generation of photon triplets using cascaded photon-pairs 
has been demonstrated in periodically poled lithium nio- 
bate (PPNL) crystals [p. This process will initiate new 
class of experiments in quantum information technologies 
using photons. Full characterization of a three-photon 
Greenberger-Horne-Zeilinger (GHZ) state using quantum 
state tomography has been performed [22j . Very recently, 
the distinction of three-photon GHZ and W states entan- 
gled in time and space has been reported by comparing 
the second-order and third-order correlation functions 
23]. As an application of this approach a method of 
generation of a narrowband three-photon W state entan- 
gled in time (or energy) via two four-wave mixing pro- 
cesses in cold atomic gas media has also been proposed 
and studied [24j . It was also shown that quantum-optical 
states based on three-photon down-conversion can be ef- 
fectively generated in cascaded optical parametric oscilla- 
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tor (OPO). The cascaded OPO was proposed [25|,[26| and 
experimentally realized by using the dual-grid method of 
quasi-phase matching (QPM) [32]. This scheme that in- 
volves cascading second-order nonlincarities is based on 
parametric processes of splitting and summing in which 
the frequencies between the pump and two subharmon- 
ics frequencies are in the ratio of 3: 2: 1. This cascaded 
configuration is different from the ones proposed in the 
cited papers. It involves the fundamental mode driven by 
an external pump field at the frequency u; an d two sub- 
harmonic modes u>2 — ^p, and oj\ = A remarkable 
feature of OPO is comparatively low generation threshold 
25j in comparison with the scheme of direct intracavity 
three-photon down-conversion, where the pump power 
threshold is determined by third-order susceptibility fl7l | . 
Another remarkable feature is formation of three phase 
locked states for both subharmonic modes equally spaced 
byf. 

In this paper, we consider semiclassical- and quantum- 
properties of photon triplets generated in cascaded para- 
metric process noted above: loo — > lj\ + L02 and L02 — > 
uji + lo\ . This setup offer a variant of the scheme realized 
in the experiment 2l| and leads to generation of triple 
photons with approximately equal frequencies. In this 
way two physical cases will be considered: (a) cascaded 
spontaneous parametric three-photon splitting in a "su- 
perlattice" of nonlinear and anisotropic linear segments 
pumped by an ultrashort pulse; (b) amplification of these 
processes in the presence of an optical cavity. 

Thus, in the first part of this paper we propose de- 
tailed studies of quasi-phase matching in cascaded con- 
figurations considering complex nonlinear structures for 
simultaneously phase matching of two parametric pro- 
cesses in optical lattices. We consider dual-grid structure 
that involve periodically poled crystals for effective gen- 
eration of photon triplets as in the experiment [211 ] . as 
well as we propose the analogical scheme for preparation 
of the three-photon GHZ polarization entangled states. 
We also consider dual-grid structures with second-order 
nonlinear and linear materials for generation of photon 
triplet with arbitrary joint spectrum. In this way we 
analyze the production of heralded two-photon polariza- 
tion entanglement in such superlattice by using the con- 
ditional method of detection of auxiliary photons. For 
this goal the various cases of the spectral factorization of 
three-photon amplitude are analyzed in details. 

The second part of the paper is devoted to studies of 
cascaded OPO on the base of the above analysis. In the 
cited papers [25|, [27j the general theory of OPO has been 
developed, however, without any concretization of the ef- 
fective coupling constants between modes subharmonics 
that are characterized by the wave- vector mismatches of 
the parametric processes ujq — > UJ1+UJ2 and u>2 — > wi+wi. 
Here, we consider these points in details. Besides this, 
the lasing of " three-photon" mode is investigated in tran- 
sition through the generation threshold. This investiga- 
tion includes also analysis of quantum statistics of modes 
i.e. photon number distributions, third-order correlation 



functions and the Wigner function for short interaction 
time as well as for over transition regime of OPO, when 
dissipative effects are included. 

The paper is arranged as follows. In Sec. II we derive 
the Hamiltonian of the system and derive three-photon 
state through QPM concurrent nonlinearities. In Sec. 
Ill we briefly describe three-wave type-II interaction in 
nonlinear superlattice structures, analyze generation of 
photon triplet in periodically poled nonlinear crystals 
(PPNC) and propose the scheme for preparation of GHZ- 
states. In Sec. IV we derive the conditions for spectral 
factorization of the three-photon amplitude. Section V 
is devoted to the analysis of cascaded OPO. We present 
the results on photon-number distributions, third-order 
correlation function of photon-numbers and the Wigner 
functions for subharmonic modes of OPO for both short 
interaction time and steady-state regime. 



II. THREE-PHOTON STATE THROUGH 
QUASI-PHASE MATCHING CONCURRENT 
NONLINEARITIES 

In this section, we consider multiple optical interac- 
tions in nonlinear media, leading to the simulta- 
neous quasi-phase matching (QPM) of two parametric 
processes. We investigate in details the cascaded pro- 
cesses of three-photon splitting and summing that lead 
particularly to the three-photon down-conversion. In 
the scheme the pump field at the frequency loq, which 
converts to the subharmonics at the central frequencies: 
cji = ^ and 0J2 = -f 1 throughout two cascaded pro- 
cesses: Wo — > oji + w 2 and u>2 0J\ + o^i- Note, that a 
complete description of photon triplet requires that all 
degrees of freedom associated with the quantized elec- 
tromagnetic field of the occupied optical modes must 
be taken into account. The specification of the photons 
involves their polarization, their wave-vector k and fre- 
quency or their spatial and time dependence. The dis- 
persion relationship between k and w reduces one de- 
gree of freedom, such that dependence from the trans- 
verse wavevectors k±_ should be taken in account. In 
this way, the transverse correlations in tripartite entan- 
glement have been studied in the paper [28| . Complete 
description in which the source produces three-photon 
states in different modes entangled in time and space has 
been done in the cited paper [23| . Analogous multimode 
calculations for cascading parametric systems in the lay- 
ered media seem to be complicated and can be carried out 
only by numerical means. However, the full discussion of 
all degrees of freedom is beyond the scope of this pa- 
per, for the further discussion we assume that the spatial 
modes of the PDC photon pairs are independently decor- 
rclatcd. We proceed from the single-mode approximation 
that selects only a single transverse spatial mode for the 
purpose of this paper. Thus, we consider collinear, one- 
dimensional configurations focusing on studies of spectral 
correlations and polarization entanglement. Such spatial 
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decorrelation can be achieved, for example, by waveg- 
uides @ , in which the photons may be emitted only into 
specific transversely-confined modes. On the other hand, 
one-dimensional collinear approach used here is also fully 
consistent for consideration of cascaded PDC, which gen- 
erates three-photon states in the presence of an optical 
cavity. This analysis will be presented in Sec.V. Thus, 
in the discussed approach the Hamiltonian of three- wave 
interaction reads as 



H i t ) = S Y, / dz X ™ fi Jz)E a (z,t)Ef l (z,t)E^z,t) = 

Ht-\t) + HW(t), (1) 

(2) 

where x a p * s tne second-order susceptibility of the 
one-dimensional medium of length L. The electric field 
operator E a (z,t) is the superposition of the three com- 
ponents corresponding to the pump field EL tCX (z,t) and 
the fields of two subharmonics consisting of positive 



Ej + J(z,t) and negative Ej~ a '(z,t) frequency electric field 
operators (J = 1,2). For simplicity, we also omit polar- 
ization states of electromagnetic fields and polarization 
indexes of nonlinear coefficient. The polarization states 
of modes will be considered in sections III and IV. In this 
case, the field operators obey 



A-) 



E { +\z,t)=i 



2^ 



Nj(u>)aj{uj) 



(2) 



where aj(u>), (j — 1,2) are subharmonic modes opera- 
tors, while classical pump field reads as 



E[ + \z,t)=i 



i(kL (uJ.z)z-Lot) 



(3) 



In these equations, Nj(w) 



is the normal- 



ized factor, where rij(u>) is the refractive index of the 
medium at the given frequency, S is the cross section 
area of the beams, kj{z,uj) — —rij(z,uj). The fields 
have frequencies centered on carried frequencies Lij whose 
bandwidths Aw a are narrow in comparison to the central 
frequencies ujq, u)\ — ^ and u 2 — 3-- 

Note, that below we consider generation of photon 
triplet due to cascading processes but not due to 
direct x interaction. In this approach the Hamilto- 
nian (flj is written trough the elementary x^ interac- 
tions that correspond to the processes: ujq — » u)\ + uj 2 
and uj 2 — > wi + lu± while three-photon down-conversion 
ujq — > u)\ + UJ\ + (J i is calculated in the standard manner 
(see, for example, |2ll. [24. 29]) in the second-order of the 
perturbation theory on x^ 2 \ ( see j formula ([TTjO . Thus, 
the Hamiltonian of the system that describes nonlinear 
X^ interactions under the rotating-wave approximation 
and for neglecting the dependencies from the transverse 
wave- vectors is calculated as 



H(t)=H 1 (t) + H 2 (t), 



(4) 



H 1 {t) = -i— J dwdw 1 dw 2 C(w,w 1 ,w 2 )Ar 1 (cj 1 )Ar 2 (w 2 )e i( " 1+ " 2 " w) *-B L (w)a + (w 1 )6 + (w 2 ) + ^.c., 
H 2 (t) = -i— ( dJzdJidw^J^J^J^Nzi^ (5) 

27T J 



Here and upwards, we use new denotations for the fields 
operators: ai = a, a 2 = b. The Hamiltonian depends on 
the Fourier spectra of the function X {z) 

C(«,o;i,« a )= [ L dz X (2) (z)e iAk ^, (6) 
Jo 

f( w ; j wi,« 3 )= / dz X {2 \z)e iAk ^ z , (7) 
Jo 

where the phase mismatch vectors are 

Aki(z) = k L (ujo,z) - ki(ui,z) - k 2 (u 2 ,z), (8) 
Ak 2 (z) = k 2 (u>2,z) - k 1 (io ll z) - ki(u)l,z), (9) 

u) = uj\ +lu 1 +lo 1 , and L is the total length of the medium. 
This Hamiltonian is presented in the general form that 
includes dispersion of interacting waves. Below it will 
be concretized for multilayered nonlinear structures with 



different susceptibilities of nonlinearity and with different 
refractive indexes, pumped by the pulse laser light beam 
with Gaussian distribution 

E l (lu) =£ exp (- ^-(lo-luo) 2 ^. (10) 

We apply now Eq.© for investigation of three-photon 
states in the cascaded processes uio w\ + u)\ The 
amplitude of this process can be calculated in the second- 
order of the perturbation theory with Hamiltonian ([5]). 
The vector state in the second-order of the perturbation 
theory is described in the Dyson series as 

m)) = (lt) j ^ j d *"tf(*')#(*")|0)a|0>6, 

(11) 

where |0) a and |0)& are vacuum states of the modes at the 
frequencies u>% and uj 2 . By using the expressions (|S|) we 
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calculate the final state of this cascaded setup for large 
time interval t — > oo. By using the commutation relations 



is calculated as 



&(cj),&+(w 



a(ui),a + (uj ) = S(u> — ui ), the result 



2h 2 



duidui^l Ni (wi)iVi (w^Ntiul , u v (J[ ) |0 b , l(w), l(w )1 (w" )) , 



(12) 



r 



where |l(k>)) = a(w)|0) is a single-mode state, 0&, l(w) are the occupation numbers of the modes and 
I 



^{w\,w 1 ,J[) = -E L (u>i +oj 1 +uj[) 



dw 2 



W 2 



(13) 



is the amplitude of three-photons at the frequencies 
u!i,uj 1 ,uj 1 , (e — ^ 0). Three-photon down-conversion is 
controlled by energy conservation between the pump and 
daughter photons cjo = wi+cjj+Wj, while the amplitude 
is determined by the pump spectral amplitude and the 
phase-matching functions. 



For the case of degenerate three-photon down- 
conversion: lo\ ~ uj 1 ~ ~ ^ the maximum of this 
amplitude is realized in the range of resonance w 2 — 
+ uj 1 — . In this case the Eq. (fT3j) approximately 
reduces to the result 



(14) 



Note, that the approximated result (fT4")l can be directly 
obtained from Eq. (fTTj) if we use the Teylor type time- 
ordering instead of the Dyson series expansion. The ba- 
sis of this approximation for the parametric processes has 
been also discussed in [3^, HH . Below, we represent appli- 
cation of these results for effective generation of photon 
triplet 21], for preparation of both three- photon GHZ 
entangled states and the heralded polarization states. 



III. THREE-PHOTON SPLITTING IN 
COMPOSITE NONLINEAR MEDIA 



In this section we first briefly discuss three-wave type- 
II interaction in one-dimensional nonlinear compos- 
ite media, consisting of alternating layers with different 
coefficients of nonlincarity. In this case, the interaction 
Hamiltonian (|T|) can be expressed as the sum of interac- 
tions in each layer in terms of the electric fields in nth 
(-)/ 



#(*) = S J2 dz X {2 \z)E; n (z,t)E[-\z,t)E! 2 -\z,t) 

n Jz " 

+h.c, (15) 

where x^ 2 \ z ) 18 the second-order nonlinearity. In ap- 
proximation that the effects of refraction from the each 
layers are small the classical pump field E pn and the pos- 
itive frequency parts of the fields of the subharmonics: 
Ej~\z,t) (j = 1,2) can be expressed as 

E pn (z,t)= f due iut E Ln {uj)e(- iko ^ z - z ^, (16) 



layer, E Ln (z,t), E) n ! {z,t) 



E)~\z,t) = J due^A* n {u)e-^^'^a+(uj). (17) 

The electric fields (US]), JT7J have been studied in HU 
for one dimensional composite materials. In the case 
of composite layered structure the coupling function: £ 
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(and/or £), (see Eqs.([6]),([7J)) could be rewritten in the 
following form (see, the paper Q) 

n Jz ™ 

m ^ ^ 

m— 1 

<Pm = 2J l nAk n , (fix = 0. (18) 
rt 

Here, Xm is the second-order susceptibility in the m-th 
layer, Z m is the length and Afc TO is the phase mismatch 
vector for the m-th layer. 

It is easy to generalize these results on the case when 
polarization structure of three- wave interaction is impor- 
tant. Note, that in this case we should not include a 
summing on the polarization indexes in the formula ([TJ, 
if the phase matching of three-waves interaction should 
take place in each mth layer only for the concrete com- 
binations of three indexes: a m , j3 m , j m . In this case, 

(2) 

the tensor of second-order susceptibility x a b 7 reduced 

(2) 

to Xm = X a 8 7 susceptibility in the mth layer. How- 
ever, in this case the phase mismatch vector for the mth 
layer Ak m also depends on the indexes a m , j3 m , j m . 
Note, that this result is convenient for realizing vari- 
ous schemes of QPM, in particular, for description of 
periodically poled nonlinear structures. The technique 
of quasi-phase-matching (QPM) and the availability of 
multiple gratings in a single periodically poled nonlinear 
crystal (PPNC) offer flexibility in achieving multiple in- 
teractions, which provide enhanced functionality of non- 
linear optical devices. 

As has been mentioned in Sec. II the maxima of the 
three-photon process depend on maxima of coupling con- 
stants of the each cascading process. Perfect QPM is the 
condition for the effective process implementation. Thus, 
in the subsection A QPM of cascading three-photon pro- 
cess: wo — > k>i + u>2, <-^2 — >• Wi + ui is considered in the 
complex nonlinear structures with dual-gratings, while 
the subsection B is devoted to investigation of frequency- 
uncorrelatcd photon triplet and subsection generation, C 
involves patterns of preparation polarization entangled 
GHZ states in such structures. 

A. QPM for generation photon triplets 

In this subsection we develop the theoretical approach 
of the cascaded generation of photon triplets, particu- 
larly, obtained experimentally in [2l|. We analyzed so 
called dual-grid nonlinear crystal in frame of mismatch- 
ing of two SPDC processes. The crystal consists of 
two periodically poled nonlinear segments with differ- 
ent periods. First segment involves M domains of the 
length 1%, with positive x ancl negative — x susceptibili- 
ties which are interchanged one to other and refractive in- 
dex rii (corresponding phase-matching function is Afci), 



and the second segment consists of N domains of the 
length I2 again with positive x and negative — x suscep- 
tibilities and refractive index 712 (corresponding phase- 
matching function is Afe). The crystal is pumped by a 
laser field, thus the first segment become primary down- 
conversion source, and photon pair is created. One of 
the photons from this pair drives the secondary down- 
conversion process in the second segment, generating a 
second pair, hence, a photon triplet. Because the photon 
triplet originates from a single pump photon, the created 
photons are strongly correlated. The ujq — > u)\ + u>2 and 
UJ2 — > u>\ + u)\ processes has been implemented in the first 
and the second sections of the crystal, correspondingly. 
We calculate the coupling constants of parametric pro- 
cesses: C and £ as a function of the finite number of non- 
linear segments. Below, we consider spectral amplitudes 
of these cascaded processes assuming the frequencies ui\ 
and ix>2 as variable quantities distributed around the cen- 
tral values u)\ = ^ and w 2 = ^p- The three-photon 

amplitude $(wi,w 1 ,w^ L ) from Eq.(fT4f is proportional to 
the product of the coupling constants and the pump field 
amplitude: El(lji +u>i), which is taken in the form 
of Gaussian distribution (fTU)) . This lattice is arranged in 
the manner that each of the Fourier functions ([6]) , (JTJ) is 
determines by one of the family of the segments. The de- 
tails of calculations based on formula (TT5|) can be found 
in [29j j . We derive the following expression for £ in terms 
of Eq. ([T8)l 



C = hxe~ 



' - . Ih A , \ sin^(Aki-qi) 

sine I — Afci ) — ; j— — ; — h 



sin l -k(Aki - qi) 



where 



sin ! f(Aki 



Aki = k L (uJo) - fci(^i) - k 2 (uJ2)- 



(20) 



Here 9i = 92 = are the harmonic grating wave 
vectors for the first and the second segments, respectively, 
d\ = 2Zi, g?2 = 2I2, L\ = Mli and L2 = NI2 are the total 
lengths of first and second sections, correspondingly. 

We assume that effective uj — > uo\ + uj 2 process goes 
in the first section, so the conditions of the perfect QPM 
is Afci = q\. Correspondingly, it is easy to see that the 
second term in (fl9| could be neglected proceeding from 
the QPM condition we are required. 

Let us write the approximation of (p~9|) for M 3> 1, 
L\ ^> l\ and Afci ~ gi. We represent phase-matching 
function in following form Afci = Ak^ + Ski, where 
Afc^ is the 0-order term of the Taylor series and Ski is 
the function of the high-order terms. We rewrite QPM 
condition in the following form Afc^ -* ~ qi, here is as- 
severated that the constant term could be compensated 
by the grating wave vector qi . Using this approach and 
assuming that liSki is very small and could be neglected 
we could find that the sine function in the (1191) becomes 
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a constant multiplier equal to — . Thus, the final approx- 
imated form of the £ is represented in following form 

C = -Mhxe-t^sinc [^{^ - q l ) S j . (21) 

This result is transformed to an expression that usually 
used for PPNC in the standard treatment. 

The Fourier function ([7]) is calculated in the following 
form 

. fh \ sin^(Ak 2 - Ql ) 

£ = hre 2 sine —Ak 2 r h 

X V2 V sin l f(Ak 2 - qi ) 



L? A k9 . . fl 2 \ sin^f-{Ak 2 -q 2 ) 
hxe 2 sine -Ak 2 , -(22) 



2 / sinf(Ak 2 -q 2 ) 



where 



Akn 



k 2 (uj 2 ) - hi^i) - h(u'i)- 



(23) 



Thus, we assume that effective u 2 



u> 1 + cj 1 process 



goes in the second section, so the conditions of the perfect 
QPM is Ak 2 = q 2 . It is easy to see that in this case 
the first term in (|2"2"j) could be neglected proceeding from 
the QPM condition we required. We could write the 
approximated form for £ similar to Eq. (|21[) . 

So we calculated the coupling constants for the both 
SPDC processes and found the QPM conditions for their 
perfect implementation. The both conditions depend on 
grating wave vectors, consequently, depend on the peri- 
ods of the segments, which are easily manageable. 



B. Frequency-uncorrelated photon triplet 

In this subsection we shortly discuss generation of 
frequency-uncorrelated triplets of photons, in analogy 
with the cascaded scheme realized experimentally in [21| . 

Thus, we assume that photon frequency-uncorrelated 
pairs in three- wave cascaded process uj — > u)\ + uj 2i 
lu 2 — > loi + loi can be generated. Note, that in typical 
SPDC experiments a monochromatic pump beam is used. 
In this situation, the sum of the frequencies of the signal 
and idler photons is fixed, thus the frequencies of photons 
are anticorrelated. A number of techniques have been 
proposed for creating frequency-correlated SPDC under 
short laser pulses This point is discussed in details 
in Sec. IV for the cascaded scheme. In the frequency- 
uncorrelated approach, we can omit the frequency de- 
pendence in the phase-mismatch functions (j6j), (J7J) con- 
sidering them on the fixed central frequencies ui, u\ and 
lu 2 with negligible band widths Au <C u>, Aui < w;. In 
this case the effective Hamiltonian ([5]) is described in the 
standard way by discreet creation and annihilation op- 
erators: a + ,b + and a, b, respectively. The commutation 



relations are [a, a + ] = 1, [6, b + ] = 1, but without using 
continuous- variable operators as in the Hamiltonian ([3]). 
For the case when photons are spatially separated the 
result reads as 



H = H 1 +H 2 , 
H x = ih(('E Q b + af - ('*E*b ai ), 
H 2 = ifr(£, a 2 a^b — £ *a 2 a 3 b + ). 



(24) 
(25) 
(26) 



Here, the subscripts label the spatial mode of photons at 
the frequencies wi , E describes the amplitude of driving 
field, C and £ are nonlinear coupling constants that cor- 
respondingly coincide with £ and £ for fixed frequencies 
uj,uji,uj 2 , accurate within normalization multipliers. We 
suppose the initial state of the cascaded system is a vac- 
uum state | Vo) = |0)&|0)i|0)2|0)3. Thus, time evolution 
of the system on the formula (fTTj) leads to the final state 
in the following form [2l| 

M*)> = -^ffi|# m = E o a\0 b , li, 1 2 1 3 ). (27) 

Here the parameters C = *C ; f — t£ describing the cou- 
pling strength between the interacting fields. 



Production of GHZ polarization entangled 
states 



Below we apply the results obtained for preparation 
of three-photon GHZ states in the cascaded scheme. 
It is possible for the case when cascaded processes in- 
volve polarized photons. Thus, we modify the above re- 
sults considering three-wave interaction in the following 
form T, a ,0 n Xa!p n ( z ) E o l (z,t)E^(z 1 t)E 7 (z,t) (see, for- 
mula ([U), where a,/3,j are the indexes of polarization 
states. 

Including into consideration the polarization states of 
the photons we assume that the type-II process ujq — > 
u>i + uj 2 create the pair of photons with a vertical V and 
horizontal H polarizations. If the pump field is oriented 
at 45° to the horizontal and vertical axes two processes 
uj — > uixiy) + ui 2 (H) and cjq — > ^i{H) + uj 2 (V) take 
place in the first crystal with effective coupling constant 
X- The second, type-I crystal is arranged in the man- 
ner that the following process: to 2 (V) — > CJ\(V) + uj 2 (V) 
and uj 2 (H) — > u->i(H) + lu 2 {H) should be realized (cou- 
pling constant k). For simplicity we will restrict our at- 
tention considering frequency-uncorrelated three-photon 
states and assume that the process under photons with 
(V) and (H) polarizations are described by the equal cou- 
pling constants. Thus, we assume that photon pairs in 



three- wave processes : ujq — > uj\ + uj 2 , uj 2 



-uj 1 have 



correlations on the polarization, but not on the spectral 
lines. In analogy with Eq s. (|2"4"1) - (|2"6")l . we model the sum 
of the corresponding parametric interactions by the fol- 
lowing effective Hamiltonian 



H = H 1 +H 2 , 
Hi = ihxEoia^b^ + ajb^ 
H 2 = ihkihiaf) 2 + b 2 {a+f) + /i.e.. 



h.c. 



(28) 
(29) 
(30) 



7 



Here, a\ and b\ are the annihilation operators at ver- 
tical polarizations, while the operators a 2 and b 2 corre- 
sponds to the horizontal polarized photons of the fre- 
quencies loi and W2, respectively. 

Time-evolution of the vector state of the system is de- 



scribed by the following formula (TTTT) with the Hamilto- 
nians (j2"51) - (j3"0)) . Choosing the initial state as a vacuum 
state \ipin) = |0) ai |0) a2 |0)b 1 |0)b 2 for all modes we derive 
the final state during time evaluation in the following 
form 



t 2 H 2 H 1 \iP) m =xkE 



|0) ai |0) Q2 |0) bl |0) b2 



(31) 



where x = xt an d k — kt are the coupling constants. 
We generalize this results by including the subscripts la- 
belled the spatial modes of photons at frequency uii (see, 
Eqs. ([25]), ([26]), ([27]). Then, the prepared three-photon 
polarization-entangled states can be written as 

\m) = X~kE Q (\V)\H)\H) + \H)\V)\V)) |0) 6l |0> h2 , 

(32) 

where the states \V) = ai\0) ai , \H) = a^|0) a2 present 
the vertical and horizontal polarization states of photons 
at the frequency u>x = Thus, we demonstrate that 
in this cascaded scheme triple photons can constitute the 
Greenberger-Horne-Zeilinger states of light. 



IV. QPM IN DUAL SECTION NONLINEAR 
CRYSTAL WITH LINEAR ELEMENTS 

In this section we investigate an improvement of above 
approach based on dual-grid nonlinear crystal for engi- 
neering controllable, factorable photon states. We con- 
sider multilayered structure of second-order nonlinear 
and linear materials. The key of this investigation is the 
idea of manipulating overall group velocity delay mis- 
matches between the various fields in structured media 
for engineering twin-photon states with special tempo- 
ral and spectral characteristics from. This approach is 
expanded here for the case of multistep parametric pro- 
cesses SPDC [2l| . We analyze the dual-grid structures to 
mismatch two processes. This method has been recently 
used in arranging of cascaded OPO 



32J. The modified 



crystal consists of two family of segments: M segments of 
the length l\, with positive x an d negative — x suscepti- 
bilities (corresponding phase-matching function is Afci), 



and the second is the structure, which consists of seg- 
ments of the length l 2 with x susceptibilities and ^ lin- 
ear optical spacers of the length I3 (corresponding two 
phase-matching functions are Ak 2 and An for the non- 
linear and the linear segments). The ujq — > uj\ + uj 2 and 
lu 2 — > itJi + loi processes have been implemented in the 
first and the second sections of the crystal, correspond- 
ingly (see, Fig. [1]). 

We calculate the coupling constants of parametric pro- 
cesses: C an d £ as a function of the finite number of 



X 



Si 



d 2 



FIG. 1: Scheme of dual section layered structure, Si is a first 
section with the length L\ , which involves nonlinear domains 
of length l\ with \+ > and X- < (with period di) and S2 
is the second section with the length L2, that involves both 
nonlinear layers of length I2 with x susceptibilities of the same 
sign and linear dispersive optical spacers of length Z3 (period 
d 2 ). 



nonlinear and linear segments. We again assume that ef- 
fective cjo — > ^1 + ^2 process goes in the first section, so 
the conditions of the perfect QPM is Aki = q (q = ^). 
Thus, the result for £ is exactly the same as the first term 
in (TU). 

The Fourier function ([7]) is calculated in the following 
form 



in(^) 



-Ak 2 

2 2 



(33) 



where d 2 — (l 2 + h) is a period of the second section, 
4> 2 = \l 2 Ak 2 + AK{1 2 + Z 3 )(^) and L 2 = f (l 2 + h) 
is the total length of the second section, where goes 
lj 2 — > uj\ + lj\ process with the following phase-matching 
functions (HI) 



Ak 2 = k 2 (uj 2 ) - ki^i) - fci(w 1 ), 
Ak 2 = k 2 (uj 2 ) - ki^) - /ti(w'i), 
AK = l 2 Ak 2 + l~ 3 An 2 , 



(34) 



where L 



, (i = 2, 3), l 2 + 1$ = 1. The condition of 



perfect QPM in this configuration is AK = 0. We rep- 
resent the phase-matching functions of nonlinear: Ak 2 



and linear: An segments as following Ak 2 = Ak 



(o) 



5k 2 



and Ak 2 — Ak 



(0) 



Sk 2 . We assume that each domain 



8 



is aligned such that zero-order constant of nonlinear and 
linear terms compensate each other l 2 Ak^ + I^Ak^ = 
0. The condition of perfect QPM requires the quantities 
5k 2 and 5k exhibit opposite signs: £2^2 = —hSn. We 
additionally require the following condition Akj ' = f- 
for zero-order therm of phase-matching function of linear 
segments. These conditions are realistic and easily satis- 
fying by managing the lengths of the nonlinear and linear 
segments, and managing the dispersive parameters of lin- 
ear segments during crystal growth. The perfect QPM of 
both processes brings to maxima of coupling constants of 
the corresponding processes which leads to excellent and 
effective three-photon down-conversion. 

A. Factorable three-photon states. Heralded 
two-photon polarization states 

Pure single photon states are probably the most fun- 
damental entities in quantum optics, and constitute the 
starting point for many optically-based quantum en- 
hanced technologies. A basic requirement for many key 
applications is the ability to generate reliably pure sin- 
gle photon wavepackets capable of high-visibility inter- 
ference. As we noted, single photon wavepackets may 
be generated from number-correlated pairs by condi- 
tional state preparation. For example, photon pairs pro- 
duced by the process of spontaneous parametric down- 
conversion (PDC) in a x^ nonlinear medium allow con- 
ditional preparation of spectrally factorized photon pairs. 
The generated pure photons remain polarization entan- 
gled, thus, the detection of one photon in the pair heralds 
the presence of the conjugate photon. 

Below we are going to analyze deduced results for the 
three-photon amplitude of down-converted light for the 
cascaded processes luq — !• u> 2 + oj\ and uj 2 — > lo 1 + uj 1 in 
the dual-section crystal (see, Fig. [T]) and find conditions 
for generation spectrally factorized photon states. Each 
of the generated modes has own polarization and we take 
into account their polarizations in the coupling constants 
£ and ^ and hence in the amplitude $(wi, uj 1 , w x ) accord- 
ing to the remark presented after formula <jTSj • Thus, we 
take into account the polarization indexes a, /?, 7 = o, e 
in the vector waves of the modes. To find the conditions 
of factorization we need to write the three-photon ampli- 
tude from Eq. (|14[) in the Gaussian form. For this reason 
we expand all phase mismatching vectors in Taylor series 
and take only the zero- and first-order terms. 

hAkt = ZiA/4 0) + T^ aVl + T 2 ^v 2 + T 3 , 7 ^ 3 , (35) 



l 2 Ak 2 = hAk^ 1 + t 23 v 2 + t a „vs, (36) 
1 3 Ak = 1 3 Ak 2 °' + q 2 ^v 2 + £>3,7^3- (37) 

Here: u = Uo, v x =U\- ^, v 2 = u)[ - ^, v z = J[ - ^f, 
v = uj 2 — The coefficients presented above have the 
following definitions: 

^ _, ( dk L (v ,f3) dki (i/i,a) i \_ 

h («Zi-<),(38) 

where i = 1,2, 3, a,fi = o, e, for the arbitrary polariza- 
tion of the laser field; 

(dk 2 {v,fi) dfci(z/ M ,7) \ 

" h \ dv U ^ d^~ l ^°) 

l 2 (u- 1 - , (39) 

for arbitrary j3 = o,e polarization of the mode with the 
central frequency in the nonlinear domains; 

f dK 2 (v,/3) dKx{v^i) \ 

Q ^ ~ h \ dv d^—^=°) 

h (vp 1 - v-i) , (40) 

where fi = 2, 3, j3, 7 = o, e, the /3 polarization of the mode 
with the central frequency in the linear domains is 
arbitrary. 

Here, Tj )Q are the temporal walkoffs between pumped 
and down-converted pulses (with the corresponding 
group velocities p and u~^ a ); i Mi7 and g Pi7 arc the 
temporal walkoffs between modes with the central fre- 
quency lo 2 and down-converted pulses with the cental fre- 
quency uj\ for nonlinear and linear domains, correspond- 
ingly (with the group velocities in nonlinear u^ 1 , u~ * and 

linear ,v~}y domains. Here Up ,Vp and u~*,t>~* 
are the group velocities of modes at frequencies uj 2 and 
cji, respectively). 

Taking into account that the pump pulse has Gaussian 
form {TOD and using the formulas (3]|, d33J) and ([58 > ([4T) 1) 
wc find the following form for the three-photon amplitude 



$ a ,p n (vi,V2,V3) = F aj p„(yi,i>2, v 3 )f a {v 1 )fp{v 2 )f 1 {vz)<& at p{vi, v 2 )$ a , 7(^1, ^3)^,7(^2, ^3), 



(41) 
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F a ,p^(yi,V2, v a ) = exp I - ^—(Ti^vi + T 2 ^v 2 + T 3 ^v 3 ) j exp ( - ^-(hjv 2 + i 3 , 7 ^ 3 + 62,^2 + £"3,7^3) j , (42) 



exp 





t% M 2 2 N 2 



Af 2 



20 



A' 2 



+ -^7T T 3,-y + -57t(*3, 7 + 03,7) 



80 




3 ' 



(43) 
(44) 
(45) 



$/3, 7 (^2, v 3 ) = exp 



®a,p{v\,V2) = exp 



$a,i{v\,v$) = exp 



M 2 + 1 



A^ 2 



2 ^ 



2 M m rr 



|(| T 2ta T 3j p + "|q-(*2,/3 + Q2,/3)(t 3n + Q 3 ,-y) 



V\V 2 



v 2 vz 



(46) 
(47) 
(48) 



r 



Let us specify the three-photon process of interest and 
write conditions for spectral factorization. We concen- 
trate on the discussion of eoo state. It means that in the 
first cascaded process: ujo ^ lo 2 +uj\ photon with central 
frequency ^ and polarization e is generated, while in the 

second process: u> 2 — > oj 1 + uj[ two photons with central 
frequencies ^ and polarizations o are generated. All 
three photons are strongly correlated, as they are gen- 
erated from one input photon. Thus, the photon with 
polarization e can be used as a heralded photon. The 
condition for generation of correlated photon pairs and 
the heralded photon is 

|$e,o,oK^3)| 2 = |$W(^)| 2 |^(^2^3)| 2 , (49) 

where the separate factor l^i^^i)! 2 = |/ e (^i)| 2 de- 
scribes one-photon probability and \$>o,,o(v2> v 3 )\ 2 = 
\fo{v 2 )<& ^ {v 2) v 3 )f (i> 3 )\ 2 describes two-photon proba- 
bility. 

Here are represented the conditions which have to be 
satisfied for generation of heralded photon 



M 2 
" + 10 
M 2 

To 



T p + -^-Ti, e T 2 ^ = 0, 

1 " - / /:;,. = 0. 



v 



(50) 



These conditions are reduced to one because as it is seen 
from Eq. ([38|) T 2 o = Tz . We found these conditions from 
the idea that the joint three-photon state in the Gaussian 
form (|4ip must be free of correlations between v\ , v 2 and 
v±, i>3 spectral components. 

Note, that each of conditions (|50l) in form coincides 
with the condition derived for the spectral factorization 
of two-photon amplitude in down-conversion [8|, |33[ . It 
is seen from Eq. (|38|) . that the explanation of this results 
reveals, that this condition may be satisfied either when 
or when u^l < uj a < u^],, (for ar- 



< u 



L0 



< u 



1 

l,e 



L.S 



bitrary polarization of the laser field /3). It means that 
the group velocity of pump pulse ought to lie between the 
group velocities of the generated pulses with different po- 
larizations. This method provides the optimal compen- 
sation for the walkoff effect [H, H3| . It should be point 
out that it is impossible for type-I down conversion, since 
in this process photons with similar polarization are pro- 
duced, consequently, with the same group velocities. We 
could manage the duration of the pump field t p and the 
length of the domain li for perfect compensation of the 
terms in the conditions. 

In the result of cascaded three-photon down con- 
version we could generate completely spectrally fac- 
torized pure tree photon states: |$ ei ,o(yi, v 2l v 3 )\ 2 = 

I^VOfl^Mfl^Va)! 2 - The conditions we have 
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to satisfy are Eqs. (l50l) complemented by the following 
condition 

M 2 N 2 

T p + -Y^ T2 - eT ^° + lo" (i2 ' + 02 '°^° + g3 ^ = °- (51) 

The last is found from the requirement that the joint 
three-photon state additionally must be free of the cor- 
relations between 1/2,^3 spectral components. We have 
observed the dual-section composite crystal with linear 
elements, which helped us to achieve perfect quasi-phase- 
matching for two simultaneous processes. The linear el- 
ements play additional positive role in the compensation 
of the terms in Eq.lfSTj). thus, in the factorization of the 
states. In such configurations, there are several param- 
eters which can be managed easily: lengths of the non- 
linear l\,l2 and linear I3 domains, number of linear and 
nonlinear domains with different lengths, the bandwidth 
of the pump pulse t p , as well as the dispersion coefficients 
of linear domains g^. The condition (|5ip is easy achiev- 
able, if the dispersion coefficients of the fields in the linear 
spacers g^ have the opposite sign to those of in the non- 
linear domains t^. Thus, we could arrange the lengths 
of nonlinear and linear domains in the second segment, 
as well as the numbers of domains such that the above 
condition (l5"Tj) is fulfilled. 

Thus, we have demonstrated the modification of QPM 
with linear compensator [!,[8| for the cascaded processes. 
It is shown that this scheme enables the engineering 
triplet of photons with desirable spectral properties. Us- 
ing the feature of polarization entanglement, it is possible 
to consider one of the photons as heralded photon and op- 
erate with two renaming states definitely knowing their 
polarizations. 

V. THREE-PHOTON CORRELATION AND 
WIGNER FUNCTIONS IN CASCADED OPO 

In the presence of the optical cavity the cascaded three- 
photon splitting displays new properties that are the sub- 
ject of this section. We assume that the cavity supports 
three modes at frequencies Uo, (*>i = ^ and L02 = =^ a . 
The mode ujq is the pump mode, steerable by an external 
coherent driving field at the same frequency too, while the 
modes uj\ and L02 are the modes of the generating subhar- 
monics. In this case, the cavity modes are described by 
discreet creation and annihilation operators a + , b + and 
a, b as in the Sec. II, III. 

The model Hamiltonian for the system, in the rotating- 
wave approximation, is given by 

H int = ih (Ee~^ at a+ - E*e iuot a ) + 

ih( (a a+b + - a+ab) + ih( (a +2 b - a 2 b+) . (52) 

Here, E describes the amplitude of the driving fields, £ 
and £ are the nonlinear coupling constants that are pro- 
portional to the functions £ and £ ©, ©, correspond- 
ingly. In difference from the Hamiltonian (|24[) this one 



includes the first term that describes excitation of the 
pump mode uq, by a classical pump field at the same 
frequency. 

The cascaded OPO is dissipative, because the modes 
suffer from losses due to partially transmission of light 
through the mirrors of the cavity. The master equation 
for the density operator p of the cavity modes in the 
Limbland form reads as 

dp 1 

-qT = TzlHintp] + ^ li{ 2a iP a t - a t a i- P a t a i)i 
' ' i=0,1.2 

(53) 

where ji are the cavity damping rates for the modes a;,, 
(t = 0,l,2). 

Below, we briefly review the semiclassical results on 
the model following the papers [25J, |27( ■ For the sake of 
simplicity let us first of all consider the case of high cavity 
losses for the fundamental mode 70 3> 71,72, assuming 
the pump depletion and, hence, its back influence on the 
modes (1) and (2) are completely neglected. Then the 
steady state solutions, for t ^> 7f 1 ,7 2 ~ 1 , is founded in 
above threshold operational regime e — > 1, where 

Eth — Tp-\/^lil2- The mean photon numbers of the 
modes are 

ni = ^|( e + 3y^T) 2 , (54) 

n2 = ^ U+v^-i J ' (55) 

and the phases are equal to 

$ 2vr , $ 2?r 
<Pi = g- + ^ 2 = "3 ~ ~3~ n ' n = ' 1 ' 2 ' ( 56 ) 

where $ is the phase of the driving field E = |_E|e 4 *. 

The stability analysis shows, that solutions (I5"4"l) . ([5"5j) 
for both subharmonic modes are stable in the range E > 
Eth- However, a bistable hysteresis-cycle behavior of the 
ni and U2, in dependence on the amplitude E, is realized 
in a small interval 1 < e < j^j. Note, that analysis of 
semiclassical solutions involving the depletion effects has 
been done in (27j . 

The second part of this section is devoted to the study 
of quantum statistical properties of high-intensity sub- 
harmonic mode u>i including calculations of the corre- 
sponding Wigner functions and correlation function. Cal- 
culations are performed on the framework of the quan- 
tum trajectory simulation method [35|. Some details of 
quantum states diffusion method (QSD), in application 
to study nonlinear processes in a cavity, can be found in 
the papers [H, H3| • 

A. Photon-number distributions and Wigner 
functions for photon triplets. 

It seems that the correlation between photons in the 
triplet, as well as quantum interference effects can be 
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FIG. 2: Photon number distributions (a, d, g), the Wigner functions (c, f, i) and its contour plots (b, e, h) for wi-mode and for 
different time intervals, t = 5 x 10~ 4 7~ 1 (a, b, c); t = 5 x 10~ 3 7 _1 (d, e, f); t = 2 x 10~ 2 7 _1 (g, h, i). The other parameters 

are: c'/7 = 200, |'/7 = 100 > 72 = 7i = 7- 



evidently displayed for short interaction time intervals 
much shorter than the characteristics relaxation time, 
t <C 7 _1 - We illustrate these effects numerically on 
the base of the master equation, however, in the regimes 
when the dissipation in the cavity is unessential and the 
dynamic of modes is almost unitary. For the cavity con- 
figuration presented the validity of such approximation 
is guaranteed by consideration of short interaction time 
1/C , l/£ <C t <C 1/71,2 provided that the coupling con- 
stants C and £ exceed the dumping rates for the modes. 

The reduced density operators for each of the modes 
are constructed from the density operator of the both 
modes p by tracing over the other mode: Pi (2) — 
Tr 2 (i)(p). The photon number distribution for u>\ mode 
is calculated as the diagonal elements P\{n) — (n\pi\n) 
on the photon- number states |n), while calculations of 
the Wigner function for u>i mode are performed by using 
its standard form in a Fock space 



Wi(p,9) =Y,PhmnW mn (p,6). (57) 

mn 

Here: p, 9 are the polar coordinates in the complex 
phase-space, which is determined by position and mo- 
mentum of the quadratures x = (a + a + )/v / 2, y = (a — 
a + )/\/2i, respectively, while the coefficients W m n(P:9) 
are the Fourier transforms of matrix elements of the 
Wigner characteristic function. 

Examples of both photon-number distribution func- 
tions and the Wigner functions for wi-mode are plotted 
in Fig. [5] for the different short time intervals. In Figs. 
[5] (a),(b),(c) we depict the results for t = 5 x 10 _4 7 _1 . 
Here, the cji-mode is in the vacuum state, therefore the 
corresponding Wigner function is Gaussian. After evo- 
lution of the system, for the time t — 5 x 10 _3 7 _1 one- 
photon state at the frequency u>i is appeared due to the 
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Tr[a +3 a 3 Pl (t)] 



(58) 



FIG. 3: Third-order correlation function versus £7 on the 
threshold E/E t h = 1, curve (1); below the threshold E/Eth = 
0.7, curve (2) and E/Eth = 0.4, curve (3). The parameters 
are: C, h = 0-2, £ h = 0.1, 72 = 71 = 7. 



process wo — > Wi + 0J2, then photon triplet is generated 
due to the cascading ujq — > + cj 2 — > 3oj! (see, photon 
distribution on the Fig. [5] (d)). The three developing 
arms on the Wigner function follow the semiclassical di- 
rections (see, formula (|56j) ) of the phase space (see, Fig. 
H^e), (f)). With increasing time intervals t = 2 x 10 _2 7 _1 
pronounced three-photon structure of the wi mode is dis- 
played on Fig. Etg) , (h), (i). As we see, in this case, the 
most probable values of photon numbers are separated by 
three photons. The Wigner function shows three phase- 
components with an interference pattern in the regions 
between them. We show the regions of quantum inter- 
ference in the contour plot (see, Fig. EJh)) depicting 
negative regions of the interference terms in the black. 
Note, that three-fold symmetry of the Wigner function 
and interference pattern have been demonstrated for the 
direct three-photon down-conversion in %( 3 ) media [381 ]. 
However, we note that the results presented here for the 
cascaded configuration are different on the details from 
the analogous calculation on the Wigner function in [38[ . 

Below, we consider photon-number distributions, as 
well as the Wigner functions of the both modes at the fre- 
quencies ui\ — ^ and 0J2 — -7^- in over-transient regime 
of OPO, when dissipative effect should be included. 



B. Photon-number correlations. 

In the paper [2l| the experimental verification of time- 
correlation between photons in triplet has been demon- 
strated. In this section, we study normalized third-order 
correlation function of the photon-numbers of mode 
u)\. This investigation is complemented by the consid- 
eration of photon number distributions of subharmonic 
modes. Thus, the correlation function for zero-delay time 
intervals is expressed as 



where rii(t) — Tr(a + api(t)) is the mean photon number. 
We calculate this quantity in transition near to the gen- 
eration threshold of OPO in the over transient regime, 
t 7j~ ,72~ • The results of numerical calculations, on 



the base of the master equation ([53|) with the Hamilto- 
nian (|52p . are shown in Fig|3]for three operation regimes 
of OPO. As we see, the curves (2) and (3) describe the 
strong correlation between three photons registered at 
the same time-intervals, that lead to three-photon super- 
bunching <?( 3 ) 3> 1. Such superbunching effect usually 
takes place for the direct three-photon down-conversion. 
We demonstrate analogous correlation for photon triplet 
plj generated in the cascaded scheme. This effect de- 
creases if the system moves to the range of generation 
threshold. At the threshold we have g^ = 1.2, for 
t-y ^> 1 (see, curve (1), Fig. 3). 



C. Transition through the generation threshold 

We continue the analysis of the sub-section A to in- 
clude dissipative effects in quantum distribution of gen- 
erated mode coi and 102- For this goal we concentrate on 
OPO operated in over transient regime, t 3> 7{" ,7^" . 
The results of the numerical calculations of the photon- 
number distributions P\{n) and P2{n) of two modes (1) 
and (2) are depicted on FigJU As we see, below thresh- 
old (FigQJa)) these distributions approximately describe 
the statistics of spontaneous radiation. It is evident that 
the locations of extremum of P(rt)-function, i.e. the loca- 
tions of the most probable values of n, may be identified 
with the semiclassical stable states (|54]) . (|55p in the limit 
of small quantum noise level (C /71 "C 1,£ /71 <C 1). 
With the increase of these ratios, due to the multiplica- 
tive character [39| of the noise in our system, the loca- 
tions of these extremum becomes slightly shifted from the 
corresponding semiclassical values n\ and ri2 (f54 |) . (|55|l . 
As shows such comparison, the Fig@|c) exactly demon- 
strates this situation. 

Another remarkable feature of cascaded OPO is forma- 
tion of phase-locked states for both subharmonic modes 
equally separated by 27r/3 on phase space. As follows 
from the semiclassical results (|54|) , (|55|) there exist three 
states of equal intensities, but with different phases. We 
present below the results of numerical calculations on 
phase-locking in quantum theory in the framework of the 
Wigner functions. In below threshold regime solutions 
of the semiclassical equations for amplitudes of subhar- 
monics are equal to zero, i.e. both modes are in the 
vacuum state, therefore the corresponding Wigner func- 
tions are almost Gaussian below threshold. The results 
at the threshold are shown in Fig. [3] A clear formation 
of phase-locked states in the transition through the gen- 
eration threshold is seen. The results above the thresh- 
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FIG. 4: Photon number distributions in transition through 
the generation threshold E/Eth = 0.7 for the mode (1) at 
the frequency Wi = (curves (1)) and mode (2) at the 

frequency u)2 = 2u;o/3 (curves (2)). The regimes are: (a); 
E/Eth = 1 (b); E/Eth = 1.4 (c) for the parameters: C /7 = 
0.2, <e'/ 7 = 0.1, 72 =71 =7- 

old are shown in Fig. [51 In this operational regime the 
Wigner functions display three-fold symmetry in phase- 
space in accordance with the semiclassical relations 
(|56|) . Note, that quantum interference pattern has been 
realized for short interaction time intervals (see, FigJ2](i)) 
is disappeared in over transient regime in the presence of 
the dissipations. 



VI. CONCLUSION 

In conclusion, we have investigated joint quantum state 
of three-photons with arbitrary spectral characteristics 
generated in optical superlattises. For this goal two cas- 
caded configurations have been considered leading to pro- 
duction of spontaneous photon triplet in cascaded PDC 
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FIG. 5: Wigner functions for the modes U2 = 2ojo/3 (a) and 
uii = ujq /3 (b) at the generation threshold for the parameters: 
C'/7 = 0-2, C'/7 = 0.1, 72 = 7i = 7 and E/E th = 1. 

and for generation of high-intensity mode through cas- 
caded three-photon splitting in optical cavity. Consid- 
ering dual-grid structure that involve periodically-poled 
crystals we have demonstrated that in cascaded type- 
II and type-I configuration triple photons constitutes 
GHZ entangled polarization states. We have shown that 
in short pulsed regime of cascaded PDC it is possible 
to control the three-photon joint spectra by using the 
method of compensation of the dispersive effects in non- 
linear segments by appropriately chosen linear disper- 
sive segment of superlattice. In the result the produc- 
tion of heralded joint states of two polarized photons 
has been demonstrated in such superlattice by using 
the conditional method of detection of auxiliary photon. 
Considering three-photon splitting in an optical cavity 
we have calculated photon-number distribution and the 
Wigner function of the mode for short interaction times. 
We have shown that in this regime, the Wigner func- 
tion displays three-fold symmetry in phase-space (phase- 
locking); the three developing arms follow the semiclas- 
sical directions of the phase space with an interference 
pattern in the regions between them. We have also 
investigated three-photon cascaded down-conversion in 
cavity in over-transient regime when dissipative effects 
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are essential. Calculating third-order correlation func- 
tion of the photon-number of sub-harmonic mode we have 
shown the strong correlation between three photons regis- 
tered at the same time-intervals that lead to three-photon 
super-bunching in below threshold regime. This effect 
decreases if the system moves to the range of generation 
threshold. We have also investigated the photon-number 
distributions and Wigner functions of sub-harmonics in 
transition through the generation threshold. 
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